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Abstract 

Using the QCD string approach the adiabatic potentials and spectra of 66-hybrid 
mesons are calculated. The results are compared to lattice studies. 

(N 

1. An assumption about existence of the hybrid mesons consisting of valence quark and 
antiquark and valence gluon was made for the first time by Okun and Vainshtein [I] just 



O 



O ' afterwards the QCD appeared. Last decade was a time of intensive studies of hybrid mesons 

both in experimental and in theoretical frameworks. 
^ ! In this report we study the system of static quark and antiquark and dynamical gluon 

joined to the formers by the infinitely thin string of background gluon field in analytic QCD 
^ ' string approach [2|, |S] based on the background perturbation theory 0j. In other words, 
we study spectrum of the vibrations of the confining string, or spectrum of the adiabatic 



- 1—1 

X 



Ph. potentials. We will pay our attention to the small and intermediate quark- antiquark separa- 
tions, which are directly related to the spectrum of heavy hybrid mesons. Using the adiabatic 
potentials, we will calculate spectra of 66-hybrid mesons. 

2. In the framework of the background perturbation theory we start from the propagator 
of the valence gluon in the background field B j^, 

Gfiv = (D 2 (B) 5^ + 2igF liV {B)Y\ (1) 

where D(B) is the covariant derivative depending on the field B; F fJLV {B) is the background 
field strength tensor, which is related to the spin effects and will be considered below. In 
the formalism of the Fock-Feynman-Shwinger representation the Green function of the hybrid 
meson with the static quark and antiquark reduces to the form (see |3] and references therein) 

G h {X,X) = j™ds J (Vz g ) Xgyg exp{-K)(W h ) B , (2) 

where Wh is the Wilson loop of the hybrid meson, consisting of the trajectories Tq, Tq and 
T g of valence quark, antiquark and gluon (see Fig. 1), and K is the kinetic energy of the 
gluon. 

In what follows we will use the area law for the Wilson loop, which is confirmed by the 
lattice QCD, see e.g. 0. The area law for hybrid mesons takes the form 

(W h ) B = exp(-a(S 1 + S 2 )), (3) 
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Figure 1: Wilson loop of hybrid meson 



where Si (S 2 ) is an area of the minimal surface bounded by the trajectores of valence quark 
(antiquark) and gluon. We parametrise the minimal surfaces introducing parameters < r < 
T and < /3i )2 < 1. It is also convenient to introduce einbein fields /i(t), vi, 2 {j, (3) and write 
the Green function of the hybrid meson in the form 

G h {X,X) = J (V 3 z g ) XgUg V^Vu 1 Vp 2 exp (- £ dr l)j , (4) 

which defines the following string Hamiltonian (U, 

H = pr-L = H + ^+ [ 1 df3 1 ^+ [ l d{3 2 ^+ [' d/h £ + t dfc^, (5) 

2 Jo lv\ Jo lv 2 Jo 2 Jo 2 



where 



2v\ Jo 2u 2 

H p2 1 x 

2(n + Ji + J 2 ) + 2A(/i + J x + J 2 ) X 

: pri ) 2 T . j \ ( P r 2 ) 2 T/ j \ 2J!J : 



+ Ji) + Ml^ + Jz) + "7' 2 2 (rir2)(pri)(pr 2 ) I , (6) 



rj r 2 rfrf 



A = (// + Jx)^ + J 2 ) - J x J 2 ^ft^, J< = / d# ftviiPi), t = l,2, 

rfr 2 Jo 

r and p are the coordinate and momentum of the valence gluon, ri j2 = r±R/2, where ±R/2 
are locations of the static quark and antiquark. Einbein fields are to be excluded from the 
Hamiltonian by the conditions 

dH 5H 

-£-=0, ^7^T = 0. 7 

One can show that einbein \x plays the role of the constituent mass of the gluon, vi{(3) the 
energy density of the background field along the string, and Ji the inertia of the string. 

3. Proceed now to the calculations of the spectrum of the string Hamiltonian at small 
and intermediate quark-antiquark separations, R < 1 fm. In this region the inertia of the 
string Ji is much smaller than the constituent mass of the gluon, Ji -C fM. Neglecting J, and 
calculating extremum over \i in the initial Hamiltonian, we arrive at the Hamiltonian of the 
potential model, 

H — \/p^ + OTl + 0T 2 . (8) 
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However, it is more convenient to calculate the eigenvalues E(n) of the Hamiltonian keeping 
the einbein /x, and then minimize them with respect to the einbein, 



This procedure is of common use in the QCD string approach, and, for the ordinary mesons, 
is justified with the accuracy better than 5% (see e.g. [2]). 

Besides the nonperturbative string part considered above, the hybrid Hamiltonian should 
contain the perturbative one V c , which is dominated by the one-gluon exchange, 

3a s 3a s a s 

V c = + — . (10) 

2n 2r 2 6R v ; 

Therefore, the Hamiltonian with the string inertia neglected takes the form 

2 

H=£ + £- + ar 1 + <rr 2 + V c . (11) 

Despite its nonrelativistic appearance, this Hamiltonian is essentially relativistic due to con- 
dition Q. 

The Hamiltonian of the string correction of the order of J/ \i derived from (0) , © has 
transparent form, 

# string = ~ (-Ll + -Ll) , (12) 
6/r Vri r 2 / 

where Lj = r$ x p. 

Spin interaction in ordinary mesons with heavy quarks of mass m have been expressed 
through the number of universal uknown potentials of order 1 / m 2 using the Wilson loop in jS] . 
In the framework of QCD string approach these potentials were calculated using the method 
of field correlators, see |2j and references therein. The role of the mass m is performed in the 
approach by the einbein //, and the results therefore are applicable to both light and heavy 
quarks. 

Spin interaction of the valence gluon with the quark and antiquark in hybrid mesons is 
generated by the strength tensor term in propagator which can be rewritten as = 
i(SB) ifc , where the spin operator S acts onto the gluon wave function ^ according to (Si^/)j = 
—ieijk^k- Performing calculations similar to ones of 0, we derive the following relations in 
Eichten-Feinberg notations, 

Vf p ) = 0, V 2 {v) =e^\ where e ^ = - — , i = l,2, (13) 

and 

y(np) = _ e(np)j yi^) = ^ where e (n P ) =arj _ ( U ) 

One can verify that the Gromes relation e + V\ — V 2 = is justified for both perturbative 
and nonperturbative potentials. The corresponding Hamiltonian consists of the following 
perturbative and nonperturbative parts [3], 
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(a) E u , n u 

(b) e+ n 9 

(c) E" n 9 , A 3 

(d) e+ n u , a u 



Table 1: Quantum numbers of adiabatic levels of Fig. 2 



JJ L S ( n P) 



o I L^S L2S 



(16) 



y ri r2 



We calculate the energy levels of the initial Hamiltonian (fTT|) variationally, using the 
Gaussian wave functions of the gluon, 



The corrections (fT2 ]l -(frB ]l amount to the decrease of the total energy about 70-150 MeV 
depending on the level. 

In Fig. 2 the resulting adiabatic potential energy levels are shown in comparison with 
the lattice results |7j. Quantum numbers of the levels in the diatomic molecule notations are 
given in Table 1. The potentials are normalized to the ground QQ-level shown in the right 
bottom corner of Fig. 2(a). We can see overall correspondence with the lattice. The curves 
tend to form three multiplets with different values of the gluon angular momentum. This is a 
confirmation that the angular moment of the gluon is a good quantum number. When R goes 
to zero, the potentials increase rapidly due to the Coulomb repulsion of quark and antiquark 
in octet state of color 577(3). Note the wrong bend of IT 9 lattice level in Fig. 2(b) and absence 
of Ug lattice level in Fig. 2(c), which clearly means that these two U g levels are not properly 
resolved on the lattice. The order of levels in our approach correspond to the lattice one, 
with the exception of Fig. 2(c). Different separations between the levels of Fig. 2(a) in our 
and lattice approaches is presumably related to the influence of the Coulomb gluon-quark 
and gluon-antiquark interaction on the gluon wave function. Variational calculations using 
the Coulomb-modified wave functions are needed to verify this assumption, which will be 
performed in future publications. 

Lowest level with the quantum numbers S~, Yl g is absent on the lattice and therefore not 
shown in Fig. 2. Comparative analysis of adiabatic potentials of other approaches can be 
found in [3]. 

4. Let us calculate spectra of masses of the heavy hybrid mesons in Born-Oppenheimer 
approximation assuming that the motion of the valence quark and antiquark is slow compared 
to the motion of the valence gluon and using the adiabatic potentials. In this case an essential 
region is the one near the minimum of the potential, where the adiabatic potential can be 
replaced by the oscillatory one. Therefore the problem reduces to the calculation of the 
spectrum of three-dimentional oscillator. 

We tune the values of parameters a s = 0.225, a = 0.18 GeV 2 , m b = 4.56 GeV to reproduce 
the experimental value of the ground state of 66-meson, = 9.4 GeV, from the Cornell 




(17) 



Mr) ~ exp (-(3 2 r 2 /2) . 
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Figure 2: Hybrid potentials with corrections included (black curves) compared to lattice ones 
(color curves) in units 1/tq = 400 MeV (t~o = 2.5 GeV -1 ). QQ-distance R is measured in 
2r ~ 1 fm. States (a)-(d) are given in Table [T] Solid, dashed and dashed-dotted curves 
correspond to A = 0, 1 and 2. Solid line at right bottom corner of FigJ2fa) represents the 
Coulomb+linear potential with a s = 0.4 and a = 0.21 GeV 2 . 
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1 = 1, l QQ = (GE) 


M bbg = 


10.83 GeV 


1 = 0, l Q Q = 1 (QE) 


M bbg = 


10.92 GeV 



Table 2: The masses of the lowest fefe-exotic states in our approach 



Lattice QCD [7] 
QCD string approach [H] 



M bb ~9=l0.% GeV 
M bb ~9=l0.M GeV 



Table 3: The masses of lowest 66-exotic states with J PC = 1 + in other works 



potential, 



V bb = 2m b -i^ + aR. 



3 R 



Then we calculate the masses of bb hybrid mesons through the equation 

M bb9 = 2m b + E osc , 



(18) 



(19) 



where E osc is the energy of 6-quarks in the adiabatic potential, approximated in the vicinity 
of its minimum by the oscillatory one. 

The lowest exotic hybrid states have quantum numbers J PC = 1~ + and may be of two 
different kinds, the gluon-excited (GE) one, with the gluon angular momentum I = 1 and 
quark-antiquark relative angular moment Iqq = 0, and quark-excited (QE) one, with I = 
and Iqq = 1. Their masses are given in Table 2. One can calculate the energy of the gluon 
excitation using the mass of the GE-state from the Table, 

AM = M b ^ E) - M T(15) = 1.43 GeV. (20) 

One can see that AM <C m& so that the adiabatic approximation for 66-hybrid mesons is 
justified. The mass of c-quark is close to the energy of gluon excitation, so that adiabatic 
approximation for cc-hybrid mesons is invalid. 

The mass of the lowest exotic state of bb- hybrid meson was calculated on the lattice j7j, 
as well as in the QCD string approach in the work [Bj. The results are shown in the Table 3. 
In lattice QCD lowest quark-excited state is absent, as it was already noted above. One can 
see from the Tables 2, 3 that the mass of the gluon-excited state computed in the adiabatic 
approximation is almost the same as ours. Note that in both cases quenched approximation 
was used. One can guess that a proper account of the see quarks will not change the result 
significantly. 

The results of the work |B] were obtained using the Hamiltonian (jllj) and the method of 
einbein field. The main difference between this and ours results is due to the large constant 
subtracted in The adiabatic Born-Oppenheimer approximation was not used in jH] and 
string- and spin-corrections to Hamiltonian were not considered. 

One can find recent comparison of existing calculations of lowest exotic hybrid meson 
masses in other models in Note that ITEP sum rules and flux-tube model predict lowest 
exotic bb hybrid mass with large error bars, which include both our numbers, and bag model 
seems to underestimate the mass. Recent calculations of GE-mode in quark-gluon model [§| 
are in good agreement with both lattice and our results. The mass of QE-mode in is less 
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than GE-one, unlike ours, see Table 2. This discrepency is presumably due to poor accuracy 
of our variational procedure for lowest adiabatic level. 

In conclusion, we have shown that the QCD string approach allows to calculate the adia- 
batic potentials and masses of heavy hybrid mesons in a good agreement with lattice QCD. 
Our further tasks are use of more elaborated variational procedure for the calculation of the 
lowest adiabatic levels, as well as accurate calculations of hybrid spectra, wave functions and 
decays both within and beyond the adiabatic approximation. 

This work has been supported by RFBR grants 00-02-17836, 00-15-96786, and INTAS 
00-00110. 
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